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_ ,.iU SUMMARY 

HO 

A representation for the indirect effect of solar radia- 
tion pressure on the semi-major axis of an artificial earth 
satellite caused by passage in and out of the earth's shadow 
is examined. The theory is then used to explain the appear- 
ance of periodicities in the observed semi-major axis of the 
orbit of Telstar 2. 
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ON THE LONG- PERIOD MOTION IN THE 


SEMI-MAJOR AXIS OF THE ORBIT OF THE 
TELSTAR 2 SATELLITE 

by 

James P. Murphy 


INTRODUCTION 

The Keplerian elements of a close artificial earth satellite undergo perturba- 
tions not only due to the earth’s oblateness, but also lunar and solar gravitational 
forces, solar radiation pressure, and atmospheric drag. 

In the elements published by the Goddard Space Flight Center for the 
Telstar 2 satellite, the effects of the oblateness have been eliminated by the 
Brouwer Satellite Theory (Reference 1). Since the satellite’s perigee height 
was about 975 kilometers, it can be considered to be relatively unaffected by 
atmospheric drag. The influence of lunar and solar gravitation and direct solar 
radiation pressure has already been analyzed for long period effects in the mo- 
tion of this satellite (Reference 3). However, this satellite undergoes indirect 
solar radiation pressure effects due to the fact that the satellite passes in and 
out of the earth's shadow. This effect is best observed in the semi-major axis, 
since this element undergoes no long-period motion from any other source. 

RADIATION PRESSURE DISTURBING FUNCTION 

The disturbing function may be written as 


R 



r cos s , 


(1) 
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where s is the geocentric angle between the sun and the satellite. Thus, 
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Also, 

r cos f = a(cos E - e) 
r sin f = a(l~e 2 ) 1/2 sin E . 

Substituting Equation (2) in Equation (1) and making use of Equations (3), the 
disturbing function R becomes 


(3) 


where 
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and 


V = 
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and where 


F = F' — 
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THEORETICAL VARIATION IN SEMI-MAJOR AXIS 

The differential equation representing the variation in the semi-major 
axis is 


But, 


da 2 BR 

dt ” na BM • 
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Thus, 
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or, 


Sa R = - 2 F [s cos E + V Vl - e 2 sin e] ** , 

where E x and E 2 are the values of the eccentric anomaly of the satellite at exit 
from and entrance into the shadow of the earth, and where 

/x - G m ( 

Thus Sa R is the perturbation after one revolution. 

SHADOW CONDITION 

In order to solve for the values of the eccentric anomaly at exit from and 
entrance into the shadow, a variation of the method suggested by Kozai in 
(Reference 2) is adopted. If the unit of distance is the radius of the earth, 
then the boundary of the shadow is expressed by 


r sin s 1 . 


(4) 


We also have 


r 

•g- cos s 


-S(cos E-e)-vVl - e 2 sin E 


(5) 


and 


— - 1 - e cos E 


( 6 ) 


If we assume S, V, and e are constants during one revolution, then the 
desired values for the eccentric anomaly may be obtained. However, the 
relations above result in a fourth degree equation in either cos E or sin E. 
Thus, it is possible to solve the quartic for, say, cos E to obtain four real roots 
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between zero and one and hence have eight possibilities for the two desired 
values of E. The number of possibilities is immediately reduced to four by 
solving the quartic for sin E also and testing whether the sum of the squares 
of a root from each equation equals unity. There can be at most four such 
pairs satisfying this condition. Two of these four roots are eliminated by 
testing whether the condition that cos s is negative is satisfied. 

The two quartics in cos E and sin E may be obtained from equations (4), 
(5), and (6). If we let x = cos E and y = sin E , the quartics are 


P 0 X 4 + p x x 3 + P 2 x 2 + P 3 X + p 4 - 0 


and 


q 0 y 4 + di y 3 + <?2 y 2 + ^3 y + % ~ 0 

where 

P 0 = 7 2 + 4> 2 

P x = 2(/3y-e4> 2 ) 

p 2 = /3 2 + lay - <p 2 (l - e 2 ) 

P 3 = 2 (a/3 + e<t> 2 ) 

P 4 = a 2 - e 2 <p 2 

% = P 0 
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and where 


a = 1 + a 2 [-1 + e 2 S 2 + V 2 (l - e 2 )J 

/? = 2a 2 e (l-S 2 ) 

7 = a 2 [S 2 - e 2 - V 2 (l-e 2 )j 

0 = - 2 S Va 2 Vl - e 2 

jj. = 1 - a 2 (l + e 2 ) (l-S 2 ) 


If in solving the above quartics any complex roots or real roots of magnitude 
greater than one appear, then these roots are discarded and the process of testing 
for permissible values of eccentric anomaly continues. It is possible that there 
are no values satisfying the conditions for shadow. This indicates that the 
satellite does not enter the shadow on that revolution, and hence the value of Sa R 
for that revolution is zero. 

APPLICATION TO THE ORBIT OF TELSTAR 2 

An examination of the mean orbital elements of Telstar 2 reveals long-period 
variations in the semi-major axis, eccentricity, inclination, argument of perigee, 
and longitude of the ascending node. The mean elements referred to are the 
"double primed" variables of Brouwer (Reference 1) in which the effects due to 
the zonal harmonics through J s have been accounted for. The periodicities 
present, in all but the first of these elements have been explained Reference (3) . 
The variation in the mean semi-major axis a" of Telstar 2 over a period of 
about 540 days is indicated in Figure (1). This variation was explained by the 
indirect effect of the solar radiation pressure. A computed perturbation for each 
revolution, Sa R , was calculated and added to the sum of the perturbations from 
each of the previous revolutions. In this way an accumulative Sa (since the 
initial epoch) was obtained. In computing the perturbations between successive 
epochs, the values of the angular variables co, fi, and A. 0 were updated after each 
revolution. The values of the other elements of the sun and satellite were up- 
dated at each epoch. Explicit relations used in updating are given in the next 
section. 
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TIME (DAYS SINCE 08/18/63 11 HRS. 55 MIN. U. T.) 
Figure 1 - Semi-major Axis Versus Time 





After computing the value of Sa since the initial epoch for each of the 77 
epochs over almost 540 days, values of the semi-major axis constant a c were 
obtained by 


a 


C 


// 

a 


Sa . 


It is interesting to note that, if the perturbations were multiplied by 4/3, the 
value of a c would have remained much more constant than it did. Such a 
quantity, a k , defined by 


a 


k 



» 


would appear only to have a slight downward trend, perhaps due to a small 
amount of atmospheric drag, over the nearly 540 days of data. 

Values of a", a c , and a k over the period studied appears in Table (1) and 
graphs of a", a c and a k versus time appear in Figure (1). 

EPHEMERIS FOR THE SUN AND SATELLITE 

The mean longitude of the sun is obtained after each revolution from 
(Reference 4): 


\ @ = 279? 69668 + 36000? 76892T + 0°.00030T 2 , 


where T is the number of Julian Centuries since January 0.5, 1900. For January 
1.0, 1960 the number of Julian days since January 0.5, 1900 is 21914.5. 

Similarly, o> and Cl of the satellite are obtained after each revolution 
from 


" - w o + _t o) 

Q = ci 0 + n(t - 1 0 ) 
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Semi-major Axis of Telstar 2 



where 


and 


CO 

fl 


3 J 2 a e , v 

“4 a 2 (i- e 2 ) 3/2 U “ 3 cos i) 

-2 a 2 (l-e 2 ) 2 C ° S 1 ' 


and where (t - t 0 ) is measured in canonical units of 806.832 seconds. At each 
epoch, new values of the Keplerian elements are also used in evaluating S, V, Sa, 
and the shadow conditions. In addition, updated values of e are obtained from 


e = 23? 452294 - 0? 013013T - 0? 000002T 2 . 


CONCLUSIONS 

In computing the perturbations, the multiplier F = F' A/m must be evalu- 
ated. To calculate the presentation area A, it was assumed that the spacecraft 
is a sphere. However, this is only an approximation to its shape. In addition, 
its surface is faceted. Since at any one time only a portion of the surface of the 
satellite is exposed to sunlight, the presentation area is only a fraction of the 
total surface area. In computing the perturbations Sa and the associated 
constant a c , it was assumed that the presentation area is 1/4 of the total area. 

If we multiply Sa by 4/3, which is equivalent to considering the presentation 
area to be 1/3 of the total area, the associated constant is a k . If we assume 
that the semi-major axis decreases slightly due to drag, then inspection of the 
curve for a k in Figure (1) indicates that perhaps a larger factor for F is in order. 

Deviations of a k from a straight line— especially near the end of the curve- 
may be due to a change in the level of solar activity. Changes in solar activity 
may result in the need to adjust the Radiation Pressure Constant F' from the 
value of -4.63xl0 -5 dynes/cm 2 which was used in computing Sa. 

Since similar periodicities occur in the semi-major axes of other satellites 
above the drag region (for example Relay 1 and Relay 2), an extended analysis of 
the motion of these satellites may result in an improved value of F' and of A 
for satellites of various shapes. 

Re-radiation from the surface of the earth and the fact that the computed 
Sa is only a first approximation to the total perturbation in semi-major axis 
may also aid in explaining the residual effects present in a c or a k . 


10 


REFERENCES 


1. Brouwer, D., "Solution of the Problem of Artificial Satellite Theory Without 
Drag," Astron. J., 64(9):378-397, November 1959. 

2. Kozai, Y., "Effects of Solar Radiation Pressure on the Motion of an Artificial 
Satellite," SAO Special Report 56, Smithsonian Institution Astrophysical 
Observatory, January 30, 1961. 

3. Murphy, J. and Felsentreger, T., "An Analysis of Lunar and Solar Effects 
on the Motion of Close Earth Satellites,” Goddard Space Flight Center 
Report X-547-65-251, June 1965 (Revised November 1965). 

4. Explanatory Supplement to the Astronomical Ephemeris and the American 
Ephemeris and Nautical Almanac , Her Majesty's Stationery Office: London 
1961. 


11 


APPENDIX A 


List of Symbols 

A Effective presentation area of the satellite 
a Semi-major axis of satellite’s orbit 
a e Mean equatorial radius of the earth (6378.388 km*) 

E Eccentric anomaly of satellite 
e Eccentricity of satellite's orbit 

F* Solar radiation pressure force constant = -4.63 x 10' 5 dynes/cm 2 * 
f True anomaly of satellite 
G Gravitational constant 

i Inclination of satellite's orbit plane to earth's equatorial plane 
J 2 Zonal harmonic coefficient in earth's gravitational potential (1.08219 x 10 ~ 3 *) 
M Mean anomaly of the satellite 
m Mass of the satellite 
m' Mass of the Earth 
n Mean motion of the satellite 
R Solar radiation pressure disturbing function 
s Cosine of the geocentric angle between the sun and the satellite 
T Number of Julian Centuries (of 36525 days) since January 0.5, 1900 
x = cos E 
y = sin E 

<V @ Mean longitude of the sun 
D Longitude of the ascending node of the satellite's orbit 
fi Mean motion of fl 

co Argument of perigee of satellite's orbit 
co Mean motion of co 
e Mean obliquity of the ecliptic 

’Currently in use at Goddard Space Flight Center. 
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